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We update Standard Big Bang Nucleosynthesis (SBBN) 
calculations on the basis of recent nuclear physics compila- 
tions (NACRE in particular), experimental and theoretical 
works. By a Monte-Carlo technique, we calculate the un- 
certainties on the light element yields (^-ffe, D, ^He and 
''Li) related to nuclear reactions. The results are compared 
to observations that are thought to be representative of the 
corresponding primordial abundances. It is found that ^Li 
could lead to more stringent constraints on the baryonic den- 
sity of the universe {Qbh'^) than deuterium, because of much 
higher observation statistics and an easier extrapolation to 
primordial values. The confrontation of SBBN results with 
^Li observations is of special interest since other indepen- 
dent approaches have also recently provided Qth^ values; i) 
the anisotropies of the Cosmic Microwave Background by the 
BOOMERANG, CBI, DASI and MAXIMA experiments and 
a) the Lyman-Q forest at high redshift. Comparison between 
these results obtained by different methods provides a test of 
their consistency and could provide a better determination of 
the baryonic density in the universe. However, the agreement 
between Q,i,h^ values deduced from SBBN calculation and ^ Li 
observation on the one hand and CMB observations on the 
other hand is only marginal. 

PACS numbers: 26.35. -l-c, 98.80.Ft 

I. INTRODUCTION 

Recently, different ways to determine the baryonic den- 
sity of the universe have been exploited. Here, we use the 
usual notations where f2f, denotes the ratio of the baryon 
density over the critical density of the universe, and 77 
is the baryon over photon ratio. They are related by 
rif,ft-^=3.65x 10^ 77 with h the Hubble constant in units 
of 100 km-s^^-Mpc^^. It is now possible to confront the 
results of these different approaches to test the validity 
of the underlying model hypothesis and hopefully obtain 
a better evaluation of this crucial cosmological parame- 
ter. Three independent methods have been used so far to 



derive the baryon density of the universe, i) the pioneer- 
ing one, (Standard) Big-Bang Nucleosynthesis (SBBN), 
based on nuclear physics in the early universe, ii) very 
recently, the study of the Cosmic Microwave Radiation 
(CMB) anisotropies and Hi) the census of H (and also 
He+) atomic lines from the Lyman-a forest at high red- 
shift. 

Regarding the uncertainties attached to each of these 
methods, those related to SBBN are probably the best 
controlled. Standard BBN depends essentially on one pa- 
rameter, T], the baryon to photon ratio since the number 
of light neutrinos is essentially known. This model, in its 
standard version, has survived for many decades show- 
ing its robustness. It success in reproducing the light ele- 
ment C^He, D, ^He and Li) primordial abundances over 
a span of 10 orders of magnitude is remarkable j^]. The 
leading uncertainties come from the observations of the 
isotopes in different astrophysical sites and the way they 
are interpreted (D in particular) to estimate the primor- 
dial abundances and the insufficient knowledge of some 
reaction rates. Lithium suffers from two drawbacks i) it 
is affected more than any other light isotope by uncer- 
tainties in the nuclear reaction rates ii) the valley shape 
in its abundance versus 77 curve leads to two possible 
T] values for a given abundance. This shape is due to 
its production modes, by ^H+'^He and ^He+'^He, re- 
spectively at low and high baryonic density. The first 
difficulty could be reduced by a better determination of 
a few key cross sections, but the second one is intrin- 
sic to the calculation. Thus to remove the degeneracy 
on the baryon density, lithium should be associated to, 
at least, one other light element, deuterium for instance. 
However, the relation between Li observations and its 
primordial abundance seems more straightforward than 
for D. 

In the case of the CMB, the n^/i^ values deduced from 
observations tend to converge but their interpretations 
are probably still model dependent. The CMB analyses 
involve many parameters, principally the various energy 
densities {fltot, ^b, ^A, respectively the total density, 
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baryonic density and the cosmological constant contri- 
bution), h, the initial fluctuation spectrum index (ris), 
the reionisation optical depth (tc) and the overall nor- 
malization. The baryon density is extracted from the 
amplitudes of the acoustic peaks in the angular power 
spectrum of the CMB anisotropics. It is important to 
note that the ratio of amplitudes between the first and 
second peaks increases with fit, in contrast with all other 
cosmological parameters. Hence, the determination of 
^Ibh^ does not suffer from the cosmic degeneracy that af- 
fects r^A and and a high precision can be expected 
. However, these values are obtained in the framework 
of inflationary models that could be altered in other cos- 
mological contexts ||] . 

The third method is based on the study of the atomic 
HI and Hell Lyman-a absorption lines observed in the 
line of sight of quasars. Quasars being the brightest ob- 
jects of the universe, they can be observed at very large 
redshift. On their line of sight, atoms both in diffuse 
or condensed structures absorb part of their radiation, 
making absorption lines apparent. It allows in particu- 
lar to study the intergalactic medium via the so-called 
Lyman-a forest. This method leads to an estimate of 
the baryon content of the Universe on large scales [^,^. 
Indeed, the evaluation of flbh^ through the study of the 
evolution of the Lyman-a forest in the redshift range 
0< z <5, though indirect because of the relatively large 
ionisation uncertainties, leads to results consistent with 
the two previous methods. 

In the following, we update the Big Bang Nucleosyn- 
thesis calculations on the basis of the recent NACRE 
compilation |6| of reaction rates supplemented by other 
recent works We performed Monte-Carlo calcula- 

tions to estimate the uncertainties on light element yields 
arising from nuclear reactions alone. Similar calculations 
have been performed recently, based on a different com- 
pilation and analysis of nuclear data (NoUett and Buries 
pof , and Cyburt et al. However, here, we put the 

emphasis on Li as its primordial abundance is more re- 
liable than that of other light isotopes {D in particular) . 
Using "^Li as the main baryometer one deduces n^h^ and 
we compare it with i) with the helium and deuterium pri- 
mordial abundances, ii) other SBBN investigations and 
Hi) independent evaluations (CMB and Lyman-a forest). 



II. OBSERVATIONAL CONSTRAINTS FROM 
THE LIGHT ELEMENTS 

Here we present the selection of astrophysical obser- 
vations that we use for the determination of the baryon 
density of the Universe from the SBBN calculation. To 
estimate primordial abundances, observations are made 
on the oldest objects that are characterized by their high 



redshift z or low metallicity [J. 

The determination of the primordial *He abundances 
is derived from observations of metal-poor, extragalac- 
tic, ionized hydrogen (HII) regions. This extraction is 
difficult to the level of precision required, due to the in- 
complete knowledge of the different atomic parameters 
involved. Olive and Skillman jl^ have studied in great 
detail the systematic uncertainties, and concluded that 
the typical errors given in previous studies are underes- 
timated by a factor of about two. The extreme values 
published [0-|l5| cover the range fom 0.231 to 0.246 (in 
mass fraction), putting little constraint on models. Con- 
sequently, in this work, ^He will not be considered as a 
discriminating indicator of the baryonic density. 

Deuterium is peculiar because, after BBN, this frag- 
ile isotope, can, in principle, only be destroyed in sub- 
sequent stellar or galactic nuclear processing. Hence the 
primordial abundance should be represented by the high- 
est observed value. It is measured essentially in three 
astrophysical sites, i) in the local interstellar medium 
(present value), ii) in the protosolar cloud (4.6 Gyr 
ago) and Hi) in remote cosmological clouds on the line 
of sight of high redshift quasars (large lookback time). 
In principle, the later sample (Hi) should be the clos- 
est representative of the primordial D value, but up to 
now, the observations lead to two ranges of D abun- 
dance values. However, very recently, it has been shown 
|l6| that the D abundance on the line of sight of the 
QSO PG1718+4607 quasar cannot be determined due to 
blending between the Lyman~a and the main hydrogen 
absorption lines, contrary to a previous study Since 
this observation was the main evidence for a high D/H 
value (^ 10""* corresponding to a low rj range; see for 
instance Ref. [p^), the very high primordial D abun- 
dance seems to have lost its support and hence only one 
range of (low) D/H values remains. However, the D 
abundance data from cosmological clouds remain scarce 
and scattered (Fig. ^, upper panel). The extreme values 
deduced from the different observations [po|-^ lead to 
the interval 1.3 x 10"^ < D/H < 4.65 x lO'^ (includ- 
ing error bars). This dispersion (amounting to a factor 
of about 3), if physical and not observational, casts a 
doubt on the direct identification of the observed val- 
ues with the primordial D abundance. Alternatively, it 
could indicate that this fragile isotope has already been 



^ Metallicity represents the abundance of metals which, in 
the astrophysical language, corresponds to all elements above 
helium. As metallicity increases in the course of galactic evo- 
lution, this is an indicator of the age of an object. Abun- 
dances of common elements like Fe (or e.g. Si), are of- 
ten taken as representative of the metallicity. The notation 
[Fe/H] = \og{Fe/H)^^^^ - \og{Fe/H)^ is often used. For 
instance, [Fe/H] = -2 corresponds to 1% of the solar (0) 
metallicity. Otherwise, D/H or ^Li/H for instance, repre- 
sent the ratio of abundances by number of atoms. 
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processed in these high redshift clouds despite their low 
metallicity |pj[] . Thus, in this perspective, averaging the 
D/H abundances measured in cosmological clouds to in- 
fer the primordial value seems somewhat inappropriate. 
In addition the observations of the absorbing cloud on the 
line of sight of QSO 0347-3818 have been analyzed using 
two methods and lead to different values (stars in 

Fig. d) . This suggest that systematic errors may still be 
important. 
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FIG. 1. Observed abundances as a function of metallic- 
ity from objects which are expected to reflect primordial 
abundances. Upper panel : observed D abundances, from 
Refs. po|j2l| , p3|j2 2[p4| (stars corresponding to two analyses of 
observations of the same cloud ^2024]). Lower panel : ob- 
served Li abundances, circles and triangles from 
Ryan et al.; squares from Bonifacio and Molaro ]3Cf|. 



On the other hand, it has been shown that there 
exists a large dispersion in the local measurements (0.5 
to 4. X 10~^). This could indicate that unknown pro- 
cesses are at work to modify the D abundance at small 
scale in our Galaxy. Thus, if it is confirmed that local 
D abundances are scattered as the result of yet unknown 
physical processes, the same thing could occur in absorb- 
ing clouds at large redshift. In addition, the lowest value 
obtained at high redshift [|| [D/H = 1.65±0.35 x lO'^) 
is uncomfortably close to both the solar system and inter- 
stellar ones (respectively around 2.1 x 10~^ and 1.5 x 10~^ 
|p6| , p7|j25| ) . It is inconsistent with even the most conser- 
vative galactic evolution models since it would require a 



negligible D destruction. Note that this lowest value |2^] 
(in parenthesis in Fig. 0) is affected by a large uncertainty 
concerning the level of the continuum and the blending of 
the relevant lines |^8| . If this questionable observation is 
put aside, a trend appears in the D/H data versus metal- 
licity [Si / H] (Fig. showing a D abundance decreasing 
when metallicity increases, as qualitatively expected from 
stellar evolution |^ . (To view the trend we do not con- 
sider the ahernative analysis Q of QSO 0347-3818 to be 
consistent with the analyses of the other observations.) 
Accordingly, the true primordial D/H value should be 
obtained from extrapolation to zero metallicity. Conse- 
quently, even though, in principle D better constrains 77 
than Li (U-shape of the curve), the value of its primor- 
dial abundance is still a matter of debate. We adopted 
the highest observed value in a cosmological cloud 
assuming that lower values are the result of subsequent 
processing. 

Compared to Z?, the determination of the Li pri- 
mordial abundance from observations leaves less room 
for interpretations. Since the discovery of a plateau 
in the lithium abundance as a function of metallicity 
(Fig. |l|, lower panel), drawn for metal poor dwarf stars 
p9| , many new observations have strengthened its exis- 
tence. The fact that the abundance does not increase 
with time (metallicity) at the surface of the oldest stars 
was interpreted as being representative of the primordial 
Li abundance [^-^. As such, these measurements 
could have been affected by two processes i) a produc- 
tion related to non thermal (spallation) nuclear reactions 
(mainly a + a) in the interstellar medium, increasing the 
amount of lithium in forming stars at a given metallic- 
ity and ii) a. depletion in the envelope of these stars. 
The contribution of the first process is small at very low 
metallicity amounting typically to less than 10% at a 
metallicity [Fe/H]^-2 |||. At [Fe/H]^-l, this contri- 
bution is more significant but remains within the disper- 
sion of the data. The second one concerns the potential 
depletion of lithium by nuclear destruction and possibly 
by diffusion and rotational mixing [ p4[ . The small scat- 
ter of the data, over three metallicity decades, on the 
one hand and the presence of the even more fragile ^Li 
isotope in a few halo stars (e.g. Ref. |^) on the other 
hand, strongly limit the amount of possible depletion. 
Hence, this effect should also be within the dispersion 
of the data. Bonifacio and Molaro have deduced 
from their large observational sample a primordial value 
: \og{Li/H) = -9.762±0.012 (statistic, la) ±0.05 (sys- 
tematic). More recently, Ryan et al. | |3l] , |36t , on the basis 
of their observations, have provided a new determina- 



tion: "^Li/H = (1.23t5!:^^) X 
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Their mean value 



and (95% confidence) limits take into account all possible 
contributions from ^ Li depletion mechanisms and bias in 
analysis. But the main difference with the earlier work 
pO| is that they have taken into account a slight rise of 
the lithium abundance due to spallation reaction leading 
to a smaller primordial abundance when extrapolated a 
zero metallicity. Accordingly, we adopt their range for 
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the primordial Li abundance. 



III. CMB AND LYMAN~a FOREST 
OBSERVATIONS 

CMB anisotropy measurements give independent es- 
timates of the baryonic density of the Universe. The 
first determinations of Qth^ from BOOMERANG and 
MAXIMA [|7|-|39| yielded r^b/i^^O.OS. The Cosmic Back- 
ground Imager (CBI), ground based, has given prelim- 
inary results in marked contrast |Q (f7h/i^w0.009) to 
BOOMERANG/MAXIMA values. The Degree Angular 
Scale Interferometer, DASI (along with its sister instru- 
ment CBI) [Q, is one of the new compact interferom- 
eters specifically built to observe the CMB. Combined 
with the large angle measurements made by COBE, it has 
been able to reveal a significant signal in the second peak 
region and has determined fli,h^ ^0 .022^'^'^^'^ (lo")- Re- 
cently, new analyses |4|@ of BOOMERANG data have 
also led to f^^/i^ = 0.022Io;o°3 (Icr). But the situation 
is not yet settled and a wealth of new data is expected 
from future ground instruments, long balloon flights and 
especially satellites (MAP, PLANCK-SURVEYOR). 

As mentioned above, the study of the baryon content 
of the intergalactic medium evolution of the Lyman-a 
forest in the redshift range < 2; < 5 leads to an evalu- 
ation of Q.hh? . Such analyses have lead to S7b/i^>0.0125 
|,|| and <0.03j|,|,||. Hui ct al. || using recent 
observations [|5|,|6|| have found f76/i2=0.03±0.01. These 
various results are summarized in Table ||, compared in 
Fig. H and discussed in Sect.^. 

IV. NUCLEAR DATA 

Most of the important reactions for Li production 
(Table are available in the NACRE compilation of 
thermonuclear reaction rates Other reactions in our 
BBN network are adapted from an earlier compilation 
[0 or more recent works In our previous studies, 

we have studied the influence of individual reactions fl^ ] 
or extreme yield limits [ p8[ obtained when considering 
all combinations of low and high rates. Here, instead, 
we have performed Monte-Carlo calculations to obtain 
statistically better defined limits, as we did in a previ- 
ous work that was limited to the NACRE reactions. 
Here we update these calculations by taking into account 
uncertainties on the remaining reactions. In the follow- 
ing, we discuss the origin of these calculated uncertainty 
limits. 

One of the main innovative features of NACRE with re- 
spect to former compilations |Q is that uncertainties are 
analyzed in detail and realistic lower and upper bounds 
for the rates are provided. Using these low and high rate 
limits, it is thus possible to calculate the effect of nuclear 
uncertainties on the light element yields. Recently, two 



other SBBN calculations have been performed: one based 
on the NoUett and Buries compilation (hereafter NB) 
and another on a partial reanalysis the NACRE data by 
Cyburt et al. Q (hereafter CFO). The NB and NACRE 
compilations differ in several aspects. The NB compila- 
tion addresses primordial nucleosynthesis while NACRE 
is a general purpose compilation. Consequently, NB con- 
tains a few more reactions of interest to BBN and a few 
more data in the energy range of interest. Also, from the 
statistical ■point of view, the rate uncertainties are better 
defined in NB, however the astrophysical S'-factors^ are 
fitted by splines which have no physical justification and 
can produce local artifacts by following to closely experi- 
mental data points. On the contrary, the NACRE compi- 
lation spans wider energy ranges, and over these ranges, 
the S'-factors are fitted to functions based on theoretical 
assumptions. 
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FIG. 2. Upper panel: astrophysical S-factor for the 
^H(p,7)^He reaction (adapted from NACRE Q, for references 
see NACRE except for SC96 [^.) Lower panel: relative dis- 
persion of residuals (Scxp/Sflt). 

Fig. H shows that outside of resonances a simple fit 
(second order polynomial in this case) is sufficient to ac- 
count for !=a2 orders of magnitude variation in S. For 
the D{p,j)'^He reaction, as in Ref. [|l8|, we use the re- 



^The (astrophysical S'-factor is defined by (t(_E) = 
exp (— 27r?7) = ^^-^ exp where here rj 

(= ZiZ2e'^/hv) is the Sommerfeld parameter. It reduces 
the strong dependency of the cross section (a) at low en- 
ergy by approximatively correcting for the penetrability of 
the coulomb barrier. 
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vised data of Schmidt et al. Ultl (SC96 in Fig. |) not 
considered in NACRE. The lower panel displays the ra- 
tio between the experimental and fitted S values. The 
dispersion around unity is small and can hardly be con- 
sidered physical. It should be noted that by using such 
a simple fit (when permitted by theory), precise data 
points outside the range of BBN energy (e.g. the high 
energy data point in this case) helps constrain the fitted 
S'-factor in the region of interest (~0.1 MeV in this case). 
This eliminates spurious local effects induced by a few er- 
ratic data points associated with experimental problems 
rather than a genuine physical effect. For instance, Fig. ^ 
(adapted from NACRE [^) shows the nuclear data and 
the NACRE recommended 5-factor for the ^ Li( p,a)'^He 
reaction can be compared with Fig. 12 of NB [Q. In this 
latter analysis, the S'-factor is unduly influenced by the 
Harmon data which is a measurement relative to the 
assumed constant S'-factor of the ^Li(p,a)^He reaction. 
The rise of S at low energy is also more likely interpreted 
by the effect of atomic electron screening of the nucleus. 
In this energy range, the cross section is expected to be 
non resonant, mainly determined by the tails of higher 
energy resonances. This is why, NACRE has constrained 
S to follow a low order polynomial with the consequence 
that the good and extensive data provided by Rolfs and 
Kavanagh are better taken into account. 



TABLE 1. Influential reactions and their sensitivity to nu- 
clear uncertainties for the production of *He, D, ^He and ^ Li 
in SBBN. 



Reaction \ AN/N 




D 


•"He 


■'Li 




n.s. 


n.s. 


n.s. 


0.08 


^H(p,7)^He 


n.s. 


-0.19 


0.19 


0.26 


^H(d,n)^He 


n.s. 


-0.09 


0.06 


0.12 


"H(d,p)^H 


n.s. 


-0.03 


-0.04 


0.01 


^H(d,n)'*He 


n.s. 


n.s. 


n.s. 


-0.07 


'H(a,7)'Li 


n.s. 


n.s. 


n.s. 


0.24 


^He(n,p)^H^^''' 


n.s. 


n.s. 


-0.06 


-0.03 


^He(d,p)''He('=''' 


n.s. 


n.s. 


-0.12 


-0.12 


''He(Q,7)'Be 


n.s. 


n.s. 


n.s. 


0.39 


'Li(p,a)*He 


n.s. 


n.s. 


n.s. 


-0.25 


^Be(n,p)^Li('=''* 


n.s. 


n.s. 


n.s. 


-0.13 



AN/N=Nh/Ni - 1 

n.s. : not significant (|AiV/iV| < 0.01). 

Chen & Savage §; Brune et al. §; Smith, 
Kawano & Malaney |^]. 
'^^ ±1ct variation. 



In some cases, [e.g. the D(d, n)'^He reaction] the NB 
compilation provides more data points in the region of 
interest than NACRE. However, considering a wider en- 
ergy range, NACRE relies on an interpolation between 
high and low energy data. It is difficult to further com- 
pare the reaction rates obtained in both compilations be- 
cause NACRE provides reaction rate limits (and in some 
cases S-factors) that can be used for e.g. subsequent 
Monte-Carlo calculations, while in NB the rate calcula- 
tion and Monte-Carlo cannot be disentangled. Indeed, 
in NB, the Monte-Carlo procedure is not applied to the 
rates but to the data points within experimental errors 
followed by spline fitting. This method is expected to 
take better into account experimental errors but is dif- 
ficult to evaluate especially because it could depend on 
the partitions of the energy interval used for spline fit- 
ting which are not given. Nevertheless the final results 
(i.e . '^He, D, ^He and ''Li yields) are in good agreement 
pof showing that both approaches are valid. CFO jl^] 
have reanalyzed the compiled data from NACRE using 
the same S'-factor energy dependences but leaving the 
scaling factors free. They found that global normaliza- 
tion factors were slightly different from the NACRE ones. 
As the fitting procedure is straightforward, the origin for 
this difference is difficult to interpret. One possibility is 
that a few data points were excluded by NACRE due to 
suspected experimental problem or physical bias (screen- 
ing at low energy). CFO subsequently determined 
scaling factors for each experiment to take into account 
systematics. From the dispersion of these factors they ob- 
tained a better evaluation of rate uncertainties including 
systematic effects. These analyses underestimate the fact 
that data on experimental cross sections are in general 
of much better quality at BBN energy or above than at 
lower energies. Indeed the cross sections for charged par- 
ticle reactions drop very rapidly at low energy (Coulomb 
barrier) making experiments more and more difficult and 
hence subject to systematic errors. In addition screening 
of the nucleus by atomic electrons is known to affect cross 
sections at low energy [Q. Hence, a scaling factor ob- 
tained from a low energy measurement is likely to be 
more affected by systematic errors that another one de- 
rived from a high energy data set, even if the quality of 
the fits is the same. This could affect the calculated CFO 
systematic error contribution. 
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FIG. 3. Astrophysical S-factor for the ^Li(p,Q)*He reac- 
tion (adapted from NACRE [^). For references other than 
HA89 and R086 H, see NACRE. 



Clearly those recent compilations or analyses 
|f7|,|6 UlO|Jll[| have all improved the determination of BBN 
rates and rate uncertainties but more progress could be 
made. This would include a better theoretical determi- 
nation of the energy dependence of 5'-factors sometimes 
just assumed to be a low order polynomial. This would 
improve the constraint given by the shape of S{E). Such 
work is under way [ p5| . 

Due to the difficulty of defining a universal statistical 
method for the wide set of reactions (each with its pe- 
culiarities) and range of temperature, the NACRE rate 
limits correspond to upper and lower bounds rather than 
standard deviations. Accordingly, we assumed a uniform 
distribution for the rates between the limits (keeping the 
mean rates equal to the recommended rates). 

Four reactions of interest to BBN are not found in 
the NACRE compilation of charged particle induced 
reactions. They are the neutron induced reactions 
n^p, iH(n,7)2H, 3He(n,p)3H, ^Be(n,p,)7Li and also 
■^He(d,p)'*He. The first reaction governs the neutron- 
proton ratio at the time of freeze-out and hence directly 
the '^He primordial abundance. The main source of un- 
certainty on this rate used to be the neutron lifetime but 
its value is now precisely (886. 7± 1.9 s) known [p6[ . So 
even though improved theoretical calculations p7| may 
introduce small corrections, this reaction is now suffi- 
ciently known for BBN calculations. The following re- 
action, ^H(n,7)^H, also relies almost exclusively on the- 
ory. A new calculation including quoted uncertainties has 
been made available recently M. For this reaction, the 
uncertainties arising from the experimental input data 
(one low-energy normalization value for the cross sec- 
tion) are expected to be much smaller than those from 



theory. The errors are given by the order of the first ne- 
glected terms in the expansion |^ . To derive the reaction 
rate and its limits, we performed numerical integrations 
using the analytical formulas for the cross section and 
its calculated uncertainties As there is no way to 
determine the statistical distribution of these theorical 
errors we adopted (as for NACRE) a uniform distribu- 
tion as for the following reaction. Following discrepancies 
and lack of documentation for the '^H(p,n)'^He reaction 
and its inverse, a new and precise measurement has re- 
cently been performed The results corroborate the 
cross section provided by the ENDF/B-VI evaluation of 
neutron data [ p8| leading to an estimated uncertainty of 
5% H. For the two remaining reactions ^He(d,p)'*He 
and T3e(n,p)^Li no new measurements are available and 
we adopt accordingly the reaction rate and uncertainties 
provided by Smith, Kawano and Malaney (In the 
more recent NB analysis the reaction rates and uncer- 
tainties are not available due to the intricate coupling of 
the fitting and Monte-Carlo methods.) They performed 
an R-matrix analysis (a standard nuclear physics method 
to tackle resonant reactions) complemented by polyno- 
mial fits to the data. Their quoted la uncertainties are 
respectively 8 and 9% and we use accordingly for these 
two reactions a gaussian distribution of errors. 
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FIG. 4. Abundances of "^He (mass fraction), D, "^He and 
^Li (by number relative to H) as a function of the baryon 
over photon ratio 77. Mean values (solid curves) and 2a limits 
(dashed curves) are obtained from Monte Carlo calculations. 
Horizontal lines represent primordial *ffe, D and ^Li abun- 
dances deduced from observations (see text). For D the dot- 
ted lines represents the range of observed values (see Fig. ^) 
while the dashed lines corresponds to the adopted value of 
Ref. !§. 



V. RESULTS 

In a first step, we complemented our previous analy- 
sis limited to NACRE reactions, by calculating the 
influence of individual reaction rates on *He, D, '^He 
and ''Li yields. Then we calculated the maximum of the 
quantity AN/N=Nhigh/Niow — 1 within the range of rjiQ 
variations for each of the 4 isotopes. Positive (resp. neg- 
ative) values correspond to higher (resp. lower) isotope 
production when the high rate limit is used instead of the 
low one (see Ref. jl^ ) . Results are displayed in Table |. 

Then we performed Monte-Carlo calculations with the 
rate distributions discussed above. For each 77 value, we 



calculated the mean value and standard deviation (cr) of 
the ^Lz yield distribution. The corresponding ±2cr limits 
are represented in Fig. 0. In Fig. is represented the 
likelihood fimctions for Li only [£^77)], D only [£^(77)] 
and for both [£^'^(77)]. The na confidence intervals are 
obtained by solving: In (£(77)) = In (>Cmax) — n^/^, for 77. 
To calculate £, we use the abundance distributions ob- 
tained by Monte-Carlo together with a normal distribu- 
tion associated with the adopted primordial abundances. 
Following the conclusions of Sect. |l[ we assumed that 
the primordial ''Li abundance is such that \ogC Li/ H) is 
normally distributed with mean -9.91 and standard de- 
viation a = 0.19/2 as given by Ryan et al. [Q. We ne- 
glected the asymmetry in the error bars (-9.9TTg ijg [^) 
by taking the largest because the smallest concerns the 
''Li lower limit that only affects the £' central dip, and 
not the 77 limits (Fig. flower panel). Unfortunately, the 
U-shape of the ''Li curve together with the Ryan et al. 
p6| values leads to a merging of the low and high 77 inter- 
vals. For comparison, we also show the likelihood func- 
tion obtained when using the Bonifacio and Molaro ]3C| ] 
older value exhibiting two 77 intervals. Their merging 
clearly originates from the new lower primordial abun- 
dance obtained when correcting for the apparent '' Li/H 
versus metallicity slope (see Fig. |l| and Sect. ||). From 
this curve, we obtain for flbh^, the range 0.006-0.016 
(95% c.L, ''Li only). 

To calculate C^, the likelihood function concerning 
D only, we adopted the highest observed value pG] 
(see Sect. O). The curve (Fig. ^ is centered on 
77=5x10^^ and is only marginally consistent with the 
C' one. The situation is even worse if the smaller 
values of the D primordial abundances are considered. 
The agreement between SBBN and the ''Li and D ob- 
served primordial abundances is impressive when con- 
sidering the orders of magnitude involved but remains 
only moderately good when trying to determine the 
ilfjh^ value. The global C^''' likelihood function pro- 
vide 17;,/i^=0.015±0.003 (2(t) which is reasonably com- 
patible with the BOOMERANG || {Qbh'^ =0.022^0° oh 
la), and DASI |l| {nbh^ ^0 m2to°ol'' 2c^) values. How- 
ever, as discussed above this is based on unsettled dis- 
crepancies on the ranges dictated separately by ''Li and 
D. Better agreement between SBBN and CMB has been 
claimed in other works ||: r2b/i2=0.020±0.002 (95% 
c.L). This result from the choice of smaller D primor- 
dial abundances, drives ri^/i^ to higher values but at the 
expense of compatibility with ''Li. Considering that the 
chemical evolution of D from BBN to present is not well 
known, but given that it can only be destroyed in this 
process, our choice of adopting the higher observed value 
seems justified. However, a better compatibility with 
the more reliably determined primordial ''Li abundance, 
would even favor a possible higher D abundance, imply- 
ing that D has already been partially destroyed in the 
cosmological cloud. Of course this would decrease the 
compatibility with CMB observations. 
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FIG. 5. Likelihood function for ^Li, D, ^ Li+D and 77 
ranges from CMB, Lyman-Q and other SBBN calculations. 
The horizontal lines represent Qth^ and rj intervals : 2a con- 
fidence limits (solid lines), twice the la limits (dash-dotted 
line) or statistically less well defined limits (dashed lines and 
arrows). Labels (a, /3, . . ., 5, 6) point to the last column of 
Table |ll were more details are given. The solid curve repre- 
sents £ for the adopted value of ^ Li primordial abundance 
while the dotted curve display, jC7 if the higher value of Boni- 
facio and Molaro |^ is adopted. The likelihood function for 
D (£^, dash dotted curve) is only marginally consistent with 
as shown by (dashed curve). and £7 have been 
normalized to £max = 1.] 



Padin et al. (2001) [M; no confidence interval given. 
Pryke et al. (2001) plj 

Riediger et al. (1998) M; estimated range of values 
from Petitjean (2001) PC 

Scott et al. (2000) @ and Hui et al. (2001) @; no 
details on statistical significance given. 



TABLE IL Comparison of Qbh^ from different methods. 
Limits are given for 2a except for values in italic (see table 
footnotes). The last column provides the labels for the Qbh^ 
ranges displayed in Fig. H. 



Method 




Fig. 5| 


SBBN+Obs (Li), this work 


0.006 - 0.016 




SBBN+Obs (Li-fD), this work 


0.015±0.003 




SBBN from Buries et al."' 


0.020±0.002 


Q 


SBBN from CFO''' He + Li 


0.006 - 0.017 


P 


SBBN from CFO''' low D 


0.017 - 0.023 


7 


CMB BOOMERANC^' 


U.U4,^_Q QQg 


1 


CMB from MAXIMA'*' 


0.0325±0.0125 


2 


CMB from CBP' 


^0.009 


3 


CMB from DASI-''' 


022+° °°^ 


4 


Lyman-Q *' 


0.0125 - 0.03 


5 


Lyman-Q 


0.03±0.01 


6 



°) Buries, NoUett & Turner (2001) |5|. 

Cyburt et al. (2001) 0. 
"'^ de Bernardis et al. (2001) |4|; 2cr interval 
approximated by twice the width of the Ict interval. 
'^^ Stompor et al. (2001) 
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VI. CONCLUSION 

Big Bang nucleosynthesis has been the subject of per- 
manent interest since it gives access to the baryon density 
which is a key cosmological parameter. Though indepen- 
dent methods are now available, the SBBN one remains 
the most reliable because i) the underlying physics is well 
known and ii) there is essentially only one free parame- 
ter contrary to other methods. It is worth pursuing the 
improvement of nuclear reaction rates and abundance de- 
termination of light elements, essentially D and "^Li. 

Our SBBN results, flbh'^ =0.006-0.016 based on "^Li 
only and 0.015±0.003 with D are good agreement with 
those from Cyburt et al. [pO[. The Qhh^ value derived by 



Buries, NoUett and Turner ||59l] {D and is in reason- 
able agreement with ours. These results are broadly con- 
sistent with the CMB ones (MAXIMA, BOOMERANG 
and DASI) and those obtained via the observation of the 
Lyman-a forest at high redshift (see also Cyburt et al. 
pOf). However, the SBBN values derived separately from 
'Li and D are only marginally compatible and when us- 
ing the more reliable indicator, lithium, the agreement 
with the CMB and Lyman~a values is also marginal. 

It is interesting to note that the SBBN ||l|,|^,|ll| re- 
sults derived from the two recent and independent reac- 
tion rate compilations (NACRE Q and NB Q) agree 
very well. Progress in the derivation of primordial abun- 
dances (D) are certainly needed, but improvement in the 
determination of nuclear reaction rates would also be of 
interest (^Li nucleosynthesis specifically). Concerning 
this last point, reanalysis of existing data constrained 
with improved theoretical input are under way [^5| . How- 
ever, it would be even more important that new experi- 
ments dedicated to precise and systematic measurements 
(e.g. in Refs. |53 ^) of the lesser known cross section (Ta- 
ble |) be undertaken. 
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